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f~^ , Abstract 

jsr , We present a simple method to obtain vacuum solutions of Ein- 

._!. ' stein's equations in parabolic coordinates starting from ones with cylin- 

f^ , drical symmetries. Furthermore, a generalization of the method to a 

"^ ' more general situation is given together with a discussion of the possi- 

^^. ble relations between our method and the Belinsky-Zakharov solitons- 
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generating solutions. 
PACS numbers: 04.20.-q, 04.20.Jb 

Introduction 

Cylindrical solutions [21 El S [51 have played an important role in the 
history of general relativity. In astrophysics, cylindrical solutions have been 
applied in the context of cosmic strings [Hj • In literature various techniques 
exist (3 El El El ^3 for generating solutions of Einstein's equations. In this 
paper we introduce a simple method for generating static and stationary 
"parabolic" solutions starting from ones with cylindrical symmetries. Our 
starting point is the stationary axially symmetric line element in the form 

m 

ds^ = r\e^\{dx^f + {dx^f) + p2d02] _ fi^dt - ujclcpf, (1) 

where x^,x^ are spatial coordinates, cp is an angular coordinate, t is a time 
coordinate, p is the radius in a cylindrical coordinate system and /, 7, io are 



functions of x^,x'^. The field equations ^H] for the Une element (^ can be 
written ^1] in the form: 

VV - jifa^ - cD„2) = , V^cD - |/„cl>„ = 0, 



T,pi + np2 c _ T,p2 - Upi d 

'Ap{pl + pl)^2 ' ^^-4p(pf + /,!) + 2' 



2pi2P2 + (pil - P22)PI ,_ 2pi2Pl - (/Oil - P22)P2 

(pI + pD ' (P? + Pi) 



^1 = ~^^2 , W2 = ^$1 (2) 

where a summation over a is implicit with a = 1,2, i.e. a;^,a;^, and 
subindices denote partial derivatives. Further, the operator V^, with V^ = 
^aa + ^da^ denotes the reduced tridimensional Laplacian (without (j)) up 
to a conformal factor: the bidimensional Laplacian is given by A = d^^^. 
Consider now parabolic coordinates, which in terms of Cartesian ones are 
given by 

X = Xpcoscj) , y = Xpsincf) , z = -(A^ — p"^). (3) 

Because we are interested in axisymmetric solutions, we will use polar cylin- 
drical coordinates p and z with 

p = Xp , z = ^{X'-p'). (4) 

The inverse of transformations @ are given by 



X = yz+^/fP + ^ , p= \J ypP^+z^ - z. (5) 

Now, if we write the field equations in parabolic coordinates, i.e. x^ = 
A, x^ = p, then for the operator V^ in these coordinates we have 

y' = dl, + dl^ + \^d^ + -^d,. (6) 

The same operator when expressed in cylindrical coordinates is 

V' = dl + dl + ia, . (7) 



By comparing the expressions © and it is easy to see that if we take 
a solution of the equations ((21) with / = /(/o),$ = $(/o), then also / = 
/(A),$ = $(A) or / = /(/i),<l> = ^{fJ,) are solutions. In this way, starting 
with a solution with cylindrical symmetries, we can obtain a parabolic one 
that is non polynomial when expressed in cylindrical coordinates. 
This paper is devoted to the discussion of the solutions so obtained together 
with an investigation of the limits of validity of the method. 
In section 1 we apply the method using as starting metric the Lewis jH] 
and the Papapetrou classes ^21 of solutions and discuss possible physical 
interpretations of these solutions. In section 2 we show that, starting with 
a static spatially homogenous solution with a G3 group of motion, a class of 
stationary solutions with a G2 group of motion can be obtained which con- 
tains as subclass the solution found in subsection 1.1. In section 3 we study 
the most general coordinate system permitted by our method. Finally, sec- 
tion 4 is devoted to a generalization of the method together with a study of 
the possible relations with the Belinsky-Zakharov (B-Z) solitons-generating 
solutions. 

1 Application of the method 

1.1 Generating solutions from Lewis ones 

For a first application we consider Lewis solutions j^ given by 

'^ PQ {P^p' - B^Q^p^-') ' 



e^^ = (f:^[^V-i?^QV-^ (8) 

where B, P, Q, and e are constants. 

For our purpose, the function ^ for the Lewis solutions must be a function 

of p. From equations ((2]) we deduce 

ie-l)2BPQz 

Thus, a necessary condition to map cylindrical Lewis solutions into parabolic 
stationary ones is e = 1. Another possibility is that B — > 00, but these 



solutions belong to the Levi Civita static class that will be discussed later. 
Taking e = 1 in ((H)) we obtain uj = b, where 6 is a real constant. By 
integrating the field equations ©, we obtain the solution 

f = aX,oj = b, e^T = \/A(A2 + /i2)3, (10) 

where a is a real positive constant. Also the function 

/ = «/" (11) 

is a solution of the system and for 7 we obtain the same expression given 
in (|lfl|) with X ^ fi. In this way we have used a subclass of Lewis solutions 
to obtain two solutions with parabolic- like symmetries. 
Now we analize the properties of (|lfl)l . First of all, metric (|10)) is Petrov 
type I and has an Abelian G2 group of motion with Killing vectors ^^ = 
dt,^'^ = d(f). Besides, it has a coordinate singularity at A = 0. In cylindrical 
coordinates it takes the form 



ds^ = ^ (^±yZ±^ [,^2 ^ ,^.] ^ ^ 1 



2^ (p2 + .2)8 «^^ + y^?Ti^ 



a^z+y^p^ + z^{dt-bd^f. (12) 

Since at p = {z axis) e^'*' = 1 for z > and e^"^ 7^ 1 for z < 0, we conclude 

that (|12|) is regular on the axis only for z > 0. 

This fact is confirmed by taking the relativistic invariants 

Tjabcd Tj Tjabcd;e Tj /~iabcd/^ (^'i\ 

-ft J^abcd 1 Ji- J^abcd;e i ^ ^abcd ' ' ' j {^'^) 

where R denotes the Riemannian tensor, C the Weyl one and ";" denotes the 
covariant derivative. The invariants are coordinate independent and when 
expressed in parabolic coordinates they are singular only for A = 0, A^+/u^ = 
0, i.e. on the z axis at 2: < 0. For example 

«.««.,, = ^5!(f^!±4 (14) 

4A3(A2+^2)5 



R^ ^ '^Rabcd;e — 17 — g ■ (15) 



and 

45a3(^/^f7AVjJ_6A^^ 

16(A2+;u2)xAf 

To study the behaviour at spatial infinity the spherical coordinates are most 
appropriate. When (|1U|) is expressed in such coordinates, it is easy to show 



that the hne element is not regular at spatial infinity (r — > oo) independently 
of the azimuthal angle 9, and at 9 = ir (p = 0, z < 0) independently of r. 
Now we study the physical interpretation of H10|) . As a first step note that it 
is possible to construct (see |15| and references therein) a "local" formulation 
of general relativity and thus to define, in this "local" frame, the analog 
quantities of the Newtonian theory. In the stationary case we can define 
a "standard" gravitational field G in a reference frame F adapted to the 
stationary spacetime (^ with a gravitational potential U given by / = e^^. 
Besides, in T we can define a time parameter T analogue to the time defined 
in a Galileian reference frame, except for the fact that T is not defined 
globally but only on the geodesic of the particle. In terms of the proper 
time T, with dr = y/—gaf3dx°'dxl^, we have 

dr = y/l^^dT (16) 

where z^* = -^,z = l — Sis the 3- velocity of a test particle with respect to 

T. 

Remember that, starting with the line element Q, a frame is adapted to 

the spacetime (0) if it is represented by coordinates x'* = x'(x*), t' = 

t'(x"), z = l — 3a = l — 4. Further, if v'^ is the 4-velocity with respect to 

r, then 7;" = j/"^. 

In particular, in this picture we can define the "relative" energy H of a 

test particle as H = -^S=e^ (.^ = 0) , where mo is the rest mass of the 

particle. For metric (jlUj) . the function H is 

mo r- mo 



'' = 71=^^ = 71^^- ^''^ 

Thus, the surfaces with constant energy are rotational parabolic. Note that 
if A ~ 0, then iiT ~ 0. However, it is in principle possible to have a particle in 
the orbit with v^ = 1 — aX which is therefore ultrarelativistic with H = m,o. 
Since surfaces with A = const, are equipotential, for the spacetime (jTU]) 
orbits exist with energy H = mo, i.e. the energy for a rest non-interacting 
particle. 

Further, when Xa < 1, the potential U = log^/J < and is thus attractive. 
When A = -, / = 1 and U = 0: in this case the particle has "local" energy 
H = /^" ^ of a free particle travelling in a Minkowskian spacetime with 



VT- 



speed u. Finally, for A > 1 we have / > 1 and thus the potential U becomes 

repulsive. 

Naturally, this does not mean that the source matter of (fTU]) is a paraboloid, 

but this is an indication that parabolic symmetries have something to do 



with the solution l\W\\ . To enforce this reasoning, we consider the coordinate 
transformation 3 which changes solution ()1U() into the static form 

27 2 

ds^ = ^ [dX^ + dfj.^] + ^d4>''^ - fdt'^ (18) 

with the same functions 7 and / given in (|1UJ) . This is done by performing 
the transformation: 

dt = Adt' + Bdcj)' , d(j) = Cdt' + Ddcj)', (19) 

where A, B, C, D are functions of the non ignorable coordinates A, ^. 
Imposing that (|18)) be equal to l|10|) . we obtain the set of equations 

^2 

' .CD = AB + CDuj'^ - ADuj - BCoj, 



f 



2 



f2 
2 I J To D 7~i, . , .2 7~i2 d21 



^2 

System (|20|) has three independent equations for four variables and thus 
admits solutions. However, transformations (|19|) are purely local, i.e. non 
integrable. This means that there exists a rotating reference frame such that 
the metric appears to be static, and the coordinates t' , (p' are admissible only 
on a geodetic. Further, note that system H2U() depends on ^: for the solution 

HlOf) this quantity is equal to %. Therefore, the angular velocity of an 
observer in the frame A, //, 0', t' with respect to which the metric appears to 
be static, depends only on /i (w is constant in the spacetime (fTU)) ). Thus, the 
angular velocity of the source, as "seen" from an observer with coordinates 
A, fj,, (f), t, has a shape given by a paraboloid of rotation. For example, with 
the ansatz dt = Adt' , d(j) = Cdt' + Ddcp', with Q = -^ and Q' = -^, thanks 
to ()2U|) we obtain 

^=a262_^2+'3'_^2^2 + ^2- (21) 

Expression 1)21(1 means that the source rotates with an angular velocity de- 
pending on /x. In fact, for Q' = 0, it follows that Q = ^i^-i^ i ■ Note that for 
// — s> 00, because in this limit dt' -^ dt and dcf)' — > d(/), formula (PT|) reduces 
to Q = Q'. The main difficulty to analize the nature of the source of pUj) 
is that this metric does not admit asymptotical Minkowskian coordinates. 
Obviously, similar arguments follow for the solution (|11() . This concludes 
our study of the physical interpretation of (|1U|) . 



1.2 Generating solutions from Papapetrou ones 

Another class of solutions that can be mapped into parabolic ones is given 
by Papapetrou jj^l ■ Starting from the metric (0) (Papapetrou gauge) , these 
solutions are characterized by 

/2 + $2 ^ 1. (22) 

It is easy to see that the most general solution belonging to the Papapetrou 
class with / = /(p), $ = $(^9) is 

1 + p^^ p^-^ + 1 

where P is a real constant. 

Once equations ((U are solved for the metric functions uj and 7, we obtain 

uj = Pz + (5 , e^^ = p^. (24) 

Note that ()23() has cylindrical symmetries only when (j) = constant, i.e. for 
planes passing through the z axis. Actually, for our method to be applicable, 
we only need solutions with / = f{p),^ = ^{p): from the field equations 
((21) it is easy to see that this implies that 7 = 7(p) and to = az + P, where 
a, P are real constants. Solution (^^ is "mapped" into 

By integrating the field equations © we get 

X — 

u;=lpp''+f3 , e^^ = ^^. (26) 

2 (A2 + p2)--i 

Also in this case we can obtain another solution by taking X ^ p. 
Generally, ((25|) admits a G2 Abelian group of motion, i.e. ^^ = d^,^'^ = 
dt- Since condition ((^ is again valid, the solution (|^3)) belongs to the 
Papapetrou class. 

With arguments similar to the ones used in subsection 1.1, it can be shown 
that the metric (|25j) has a coordinate singularity at A = 0, is regular on the 
z-axis only at z > and is singular for 2; < 0, i.e. in the limit /x ^ it follows 
that e^i' -^ 1. In fact, the invariants are singular for A = 0, A^ + p^ = 0. 
Solution (|25|) is Petrov type O (flat) for P = 0, is Petrov type D for P = ±2 



and otherwise is Petrov type I. While at spatial infinity (r — > oo) it is not 
regular. Besides, for the "relative" energy H oi a test particle of inertial 
mass rriQ and velocity v we get 

H = ^^e^= ^"^° ^ ^ (27) 

and thus, orbits with constant energy in a frame T adapted to the stationary 
metric (^ are again rotational parabolic surfaces. Then, for A ~ 0, i/ ~ 
^m=^/2Al and for A ^ cx), // ~ -^^^ : in both cases if ^ 0, i.e. 

Vl — I' Vl — i^ _\-2" 

U -^ — oo. Moreover, for spacetime H25|) . C/ < with equality only at A = 1: 
therefore, at A = 1 ii = ^m=. 

Obviously, also in this case, this does not mean that the source matter of 
H25() is a paraboloid, because the shape of the configuration depends on the 
"match" between the gravitational and the centrifugal force, which is not 
given a priori. In practice, the fact that orbits with parabolic symmetries 
are allowed does not guarantee that the source is a paraboloid but is a sign 
in this direction as well as the fact that the invariants of H26() are singular on 
the axis at A = 0. Finally, since for fl^ ^ = F{X, fi) , the arguments that 
lead to (|21|) are not valid, and consequently the shape of the angular velocity 
of the source, as "seen" from an observer at rest in a general spacetime point 
with coordinates A, fi, (p, t, (since the metric is not asymptotically flat, there 
does not exist a "privileged" Minkowskian observer at spatial infinity), has 
not parabolic symmetries. 

2 Stationary solutions from static ones 

To apply our method in the above section we have taken as "starting" metric 
a subclass of Lewis solutions with a G3 group of motion and a G2 subclass 
of Papapetrou solutions (a / 0). All the generating solutions have a G2 
Abelian group of motion. In this section we show that it is possible to 
obtain a stationary G2 solution starting with a static G3 solution with Killing 
vectors ^^ = dt^S,"^ = 5^,^^ = dz- As an example we start with the most 
general cylindrically symmetrical static solution with an Abelian G3 group of 
motion (found by Levi-Civita fl]) that can be obtained from Lewis solutions 
by setting B = (uj = 0) . With the same technique used above we obtain 

f = aX^ , e2-^ = ^^^^, (28) 

(A2+^2)--l 



and the solution with A — > ;U. 

When expressed in terms of p, z, with the help of ©, all these parabolic 
solutions have complicated expressions. Note that, thanks to the analoguey 
between the Laplacians © and we can "map" all static cylindrical so- 
lutions into parabolic ones. Solution (|28|) is Petrov type O for e = 0, 2 and 
otherwise is Petrov type I. Furthermore 

o A~2~ 

/ = aX^ , e "^ = 2 , CO = Q = const. (29) 

(A2+^2){x-l) 

is a stationary solution. This solution has similar features to (|1()|) . i.e. it has 
a physical coordinates independent singularity on the z-axis at z < and 
is not regular at spatial infinity. Note that for e = 1 the class of solutions 
(Enj reduces to JTUl). Further, for the solution (gH), ^ = ^A^^-s^) and thus 
the angular velocity of the source, in general, has not parabolic symmetries. 
Now, the "local" energy H in the reference F is 

and consequently —oo<U<+oo. The arguments used for (|17|) are still 

valid. 

We now consider the interesting case of the static subclass of (|28j) with e = 2. 

This is a flat solution except on the z axis at z < and at spatial infinity 

where the metric is not regular. In parabolic coordinates we have 

ds"^ = dX^ + dfi^ + fi'^d(p'^ -X'^dt'^. (31) 

Performing the spatial coordinate transformation 

X = n cos (j) , y = fj. sin (p , z = X, (32) 

the metric (|^T|) becomes 

ds^ = dx'^ + df + dz'^ -z'^dt^. (33) 

The solution admits a G4 group of motion and is isomorphic to the the flat 
static Das |16j solution, but with i > 0, and therefore covers the Das solution 
only for z > 0. In parabolic coordinates the four Killing vector are ^^ = 



(O,cos0,-^,O),e2^(O,sin<^,2^,O),e3 = (O,O,l,O),e^ = (0,0,0,1). 
Also in this case we can generalize solution ()33() to a stationary one obtaining 



ds'^ = dX^ + (i^2 + ^2(i02 _ A2((it - Qd(t)Y (34) 

where Q is a constant. 



3 Validity of the method 

In this section we characterize the most general coordinate transformation 
for which our method works. First of all note that equations © have been 
written starting with the line element (0) (Papapetrou gauge). It is easy 
to see ^1] that the Papapetrou gauge is preserved if and only if analytical 
coordinate transformations are considered, i.e. x^ + ix"^ = F{x^ + zx^). 
Obviously, the imposition of the Papapetrou gauge does not represent a 
loss of generality. Therefore we can restrict our consideration to analytical 
coordinate transformations (note that if such coordinate transformations are 
not used, equations (HJ assume a very complicated expression). 
Since the starting point of our method is the analoguey between Laplacians 
in cylindrical and parabolic coordinates, it is natural to ask if there exist 
some other coordinates u, v, such that the (reduced) Laplacian takes the 
form 

V' = dl^ + dl + ^du + ld,. (35) 

As a first step we consider separable analytical coordinate transformations: 
p = U{x^)V{x'^) , z = z{x^,x'^). (36) 



have U^iV = z^2, UV^i = —z^i: these lead to 



Since the analyticity condition p^i = z^2, p^2 = —z^i must be imposed, we 

these lead to 

V^2 fudx^ + h{x^), 



f/xi.l ^.2x2 



U V 

System H37() only admits solutions given by 



(37) 



U = sinh X , V = sin x , z = — cosh x cos x , 


(38) 


U = coshx^ , V = cosx^ , z = sinh x^ sin x^. 


(39) 


U = e^ , V = sinx'^ , z = -e^ cosx^. 


(40) 


U = x\V = x\z = l[{x'f-{x')\ 


(41) 



Solutions (|38() and H39|l represent respectively spheroidal prolate and oblate 
coordinates. While H4U|) represents spherical coordinates and ()41() parabolic 
coordinates. Coordinates (|38j) - (|41)) are the only analytical coordinates sep- 
arable with respect to the operator V^ (V^ = S^i^i + ^^^23,2 + P{x^)dx^ + 



10 



G{x'^)d^2). Note that the parabohc coordinates are the only ones with re- 
spect to which the "similarity" reasonings with the cylindrical polar coor- 
dinates are available. However, the "similarity" reasoning with respect to 
the operator V^ can be done for another pair of coordinates. In fact, if 
we consider spheroidal prolate coordinates /i, with p = sinh/xsin0, 2; = 
cosh ^ cos 0, the operator V^ is 



72 _ ^2 , ^2 , cosh^^ ^ cos_^ ^ 
sinh ^ ^ sin ( 



V^ = 5^. + 5|. + ^S^5. + ^5,. (42) 



Now, if we take spherical coordinates compatible ^1] with Papapetrou gauge 
(P3), i.e. p = e" sin ■!?, z = e^ cos -Q, we have 

y^ = dl + dl, + d, + ^^d^. (43) 

smw 

Thus, if we have spherical solutions with / = f {'&)-, $ = ^{'&)^ we can obtain 
prolate spheroidal ones with 1? — > 0. For example, we have the stationary 
solution in spherical coordinates belonging to the Papapetrou class given by 

s\T?-d 2cosi9 

1 + cos^i? ' 1 + cos^i? 

From solution (|44() we can obtain another one with ■!? ^ 0. 
What happens if we consider coordinate transformations of the form p = 
G{x^,x'^), z = H{x^,x'^) ? It is easy to show that the analiticity condition 
G^i = H^2 , G^2 = —H^i leads inevitably to the condition G{x^,x'^) = 
U{x^)V{x'^), that has been analysed above. 

As a final consideration note that, since / = /(A) and ^ = ^(A), the term 
-dfj_ in the operator V^ as expressed by parabolic coordinates acts trivially 
on /(A), *I>(A). This means that if we consider the coordinate transformation 
given hy p = XF{p) , z = |A^ ~ '2a ^ ^ ~ const., the method is again 
applicable, with the inverse given by 



A 



-i^/^T^ + z) , F=Ja{^/^T^-z). (45) 



Obviously, the Papapetrou gauge breaks down, but because of the indepen- 
dence of / and $ on /i the field equations for these functions are again of 
the form given by the first two equations of © with V^ = d'^)^ + jdx. 
However, thanks to the second of equations H45|) . it is easy to see that the 
function F represents the possibility of performing a general coordinate 
transformation for p that does not lead to a difi'erent solution. No other 
possibilities are allowed. 



11 



4 Further improvements and final remarks 

Our starting point has been the hne element (^ written in the Papapetrou 
gauge. However, the most general form of the metric for a spacetime ad- 
mitting a two-dimensional Abelian group of isometries with Killing vectors 
dt,d^ is UH] 

ds^ = he^^{{dx'^f + {dx^f) + W^dcf)'^] - f{dt - Lod^f (46) 

where W = W{x^,x'^) with W > 0. It is easy to see that the field 
equations imply 

AW = W^i^i + W^2^2 = 0. (47) 

Condition (|47|) means that W can be choosen as a coordinate. Further, note 
that the determinant of the 2-metric g spanned by the Killing vectors dt , d^ 
is 

{-det\\g\\)'^ =W. (48) 

The function W characterizes a measure of the area of the orbits jl7| of 
the isometry group and thus it has a geometrical significance. The case 
with W = x^ = p has been analysed above. Therefore, we can start with 
a given W = W{p, z) = x^ satisfying condition (|47j) and with the other 
coordinate U given by Wp = Uz , Wz = —Up and then apply our method. 
The field equations for this new gauge are again given by provided that 
p is substituted with x^. The similarity reasoning with the coordinates W, U 
is possible only with the coordinates a, (3 given by 

W = a(5 , [/= -(a2-/?2). (49) 

For example, we can take W = \ with A = afi , p = 2(0^^ ~ Z^^) ^-^d use 
as a starting metric the solution (jSJ with e = 1 and with p substituted with 
A. Therefore, by changing W{p, z), we can increase the class of coordinates 
with respect to which our method is applicable. 

As a final consideration, we analyze the possible relations between the 
method presented in this paper and the well known jlUUllj Belinsky-Zakharov 
solitons-generating technique. This method derives from the application of 
the "quantum inverse scattering method" (QISM) to spacetimes admitting 
a two-dimensional Abelian group of isometries acting orthogonally transi- 
tively on two-dimensional spacelike or timelike orbits. For the spacelike case 
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we have two spacelike Killing vectors and the Ernst equation becomes hy- 
perbolic. In what follows we restrict our attention to the timelike case. 
The starting point of the B-Z method is the line element 

ds"^ = F{dp^ + dz^) + QABdx^dx^ , A,B = 3,4. (50) 

For a given initial solution (^0,90) of the Einstein equations, by introducing 

a 2 X 2 complex matrix ip, the field equations for FqjQq can be reduced 

to a pair (Lax pair) of two first order differential equations which are the 

1 
integrability condition for the Ernst equation with {—det\ l^ol I) ^ = a , Aq = 

0. By chosing a as one of the coordinates (the determinant of the 2-metric 

spanned by Killing vectors), as the other spatial coordinate /3 we can take 

its coniugate: a^i = (3^2 , a^^ = —Px^- The B-Z method can generate new 

solutions F, g by introducing a spectral parameter (^ together with the initial 

condition -(/'(C = ^) = 9o (for more details see ^l^HEl^])- Further, 

soliton solutions can be added depending on the pole singularities of the 

matrix ip{x^,x'^,0 in the complex plane of the spectral parameter C. The 

poles of the ip matrix are a certain function //^ = /i(p, z) of the coordinates 

and the integer index k denotes the number of solitons added to the initial 

solution. In B-Z the pole trajectory, contrary to the original QISM, are not 

constant and are given by 

2 2 - 

Hk = i^k - z ± [{ujk - z) +p]^ (51) 

where the constant w^ are the origin of the z coordinates. 
Finally, for the new solution F, g we have again 

{-det\\g\\)'2 = a. (52) 

There are formal analogies between the method presented in this paper and 
the B-Z method. First of all, both of them need a starting metric and the 
group of isometrics must be at least a G2 Abelian group. Moreover, the 
area of the orbits on the isometries group (condition ()52|)). is preserved. 
Our method requires harmonic coordinates in order to achieve the conserva- 
tion in form of the line element and therefore of the field equations. For the 
B-Z method the existence of such coordinates a, /3 is a consequence of the 
integrability condition of the Ernst equation (see P^J). Moreover, in such 
coordinates the field equation in B-Z involving i/j assume the most simple 
form. However, these formal analogies are not sufficient to etablish a direct 
link between the two methods. In fact, in spite of the similarity between the 
pole trajectory formula and expression 0, we have not at our disposal a 

13 



simple way to relate a transformed solution in our method with the number 
of solitons that can be added in the B-Z method. 

However, is interesting to note that in both cases, starting with a spatially 
homogenous metric with a G3 group of motion, the resulting solution has a 
G2 group of motion ^Sj. We conclude saying that our method is also appli- 
cable for spacetimes admitting two spacelike Killing vectors. In this case the 
equation satisfied by W is the wave equation instead of (|47|) . Furthermore, 
for the coordinates W, U, the analytical condition is substituted with 

VF,i = U,2 , W,2 = U,i , (53) 

implying that W^i^i — W^2^2 = 0. 
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